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Abstract 

The zeros of the eigenfunctions of self-adjoint Sturm-Liouville eigenvalue 
problems interlace. For these problems interlacing is crucial for completeness. 
For the complex Sturm-Liouville problem associated with the Schrodinger 
equation for a non-Hermitian 'PT-symmetric Hamiltonian, completeness and 
interlacing of zeros have never been examined. This paper reports a numer- 
ical study of the Sturm-Liouville problems for three complex potentials, the 
large- JV limit of a —{ix) N potential, a quasi-exactly-solvable — x A potential, 
and an ix 3 potential. In all cases the complex zeros of the eigenfunctions 
exhibit a similar pattern of interlacing and it is conjectured that this pattern 
is universal. Understanding this pattern could provide insight into whether 
the eigenfunctions of complex Sturm-Liouville problems form a complete set. 
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I. INTRODUCTION 



The spectra of many classes of non-Hermitian PT-symmetric Hamiltonians are real and 
positive It is believed that the positivity of the spectra is a consequence of VT 



symmetry. Examples of heavily studied PT-symmetric Hamiltonians are |l]|| 

H = p 2 -{ix) N (N>2) (1.1) 

and § 

H = p 2 - x A + 2iax 3 + (a 2 - 2b)x 2 + 2i(ab - J)x (J integer, a 2 + 4b > K ciitical ), (1.2) 

where identical grows with increasing J. For Hamiltonians like those in (|1.1|) and (|1.2|) the 
Schrodinger equations for the fcth eigenfunction, 

H* k (x) = E k ^ k (x), (1.3) 

involve a complex potential and may require x to be complex for the boundary conditions 



to be defined properly |I4fl . Thus, the Schrodinger eigenvalue problems may be regarded as 
analytic extensions of Sturm-Liouville problems into the complex plane. 

The eigenf unctions of a conventional self-adjoint Sturm-Liouville problem are complete. 
Completeness is the statement that a given function can be represented as a linear super- 
position of the eigenfunctions: 

f(x) =5>A(z). (1.4) 

n 

It is necessary for the zeros of the eigenfunctions in the complete set to become dense on the 
interval in which the Sturm-Liouville problem is defined. If the zeros did not become dense, it 



would be impossible to represent a rapidly varying function |T5|fl6[] . For conventional Sturm 



Liouville problems one can prove that the zeros of successive eigenfunctions interlace, and 



this interlacing of the zeros ensures that the zeros become dense |To| . 

A major open mathematical question for PT-symmetric Hamiltonians is whether the 
eigenfunctions form a complete set. If the zeros of the eigenfunctions of complex Sturm- 
Liouville eigenvalue problems exhibit the property of interlacing, this provides heuristic 



evidence that the eigenfunctions might be complete. A proof of completeness would require 
that we identify the space in which they are complete, and we do not yet know how to do 
this. Nevertheless, if we can understand the distribution of the zeros of the eigenfunctions, 
we gain some insight into the question of completeness for eigenfunctions of PT-symmetric 
Sturm-Liouville problems. 

II. SOME EXAMPLES OF DISTRIBUTIONS OF ZEROS 

We have studied three different complex "PT-symmetric Hamiltonians. We find that in 
every case the qualitative features of the distribution in the complex plane of the zeros of 
the eigenfunctions are very similar: We observe a shifted interlacing of zeros. We believe 
that this pattern of zeros is universal. 

Example 1: Large- A limit of the —(ix) N potential. The large- A limit of a — (ix) N 
potential is exactly solvable ||. In Fig. [I] the zeros of the 14th and 15th eigenfunctions 
are plotted and clearly exhibit a form of interlacing in the complex plane. For convenience, 
we have scaled the zeros by dividing by the magnitude of the turning points; we have then 
performed a linear transformation to fix the turning points at ±1. The appropriate scaling 
is 

z = (xE- 1 ^ + i)N/it. (2.1) 

As shown in Fig. |2|, the zeros of the first 15 eigenfunctions interlace and appear to become 
dense in a narrow region surrounding an arch-shaped contour in the complex plane. This 
contour is the Stokes' line that joins the turning points; that is, it is the path along which the 
phase in the WKB quantization condition is purely real (and thus the quantum-mechanical 
wave function is purely oscillatory) (It is interesting that this WKB path differs from 
the path that a classical particle follows in the complex plane as it oscillates between the 
turning points. The path that a classical particle follows is an inverted arch-shaped contour 
between the same two turning points. 0) 
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Example 2: Quasi-exactly-solvable —x A potential. Next, consider the quasi-exactly- 
solvable potential in ( |1.2| ) with a = 10 and 6 = 2. The eigenf unctions of the Schrodinger 
equation have the form of an exponential multiplied by a polynomial. The zeros of this 
polynomial are easy to calculate numerically. For a given J the polynomials in the eigen- 
functions all have the same degree and, as a result, all of the eigenfunctions have the same 
number of zeros. However, for the kth wave function, J — k zeros lie along the branch cut 
on the positive imaginary axis, and we consider these zeros to be irrelevant. For values of J 
ranging from 1 to 21, the qualitative behavior is always the same. In Fig. || the results for 
J = 21 are plotted and the relevant zeros again lie along the WKB paths in the complex 
plane. In this case the zeros are not contained in as narrow a region of the complex plane 
as for the —(ix) potential because the zeros have not been scaled as in Fig. |2|. The zeros 
have an imaginary part that becomes more negative as k increases and exhibit the complex 
version of interlacing. 

In Fig. |] the zeros are scaled so that z = x/\xtp\, where \xtp\ are the magnitudes of 
the classical turning points. (The classical turning points are the roots of — x 4 + 20ix 3 + 
96x 2 — 2ix = E.) The scaled zeros lie in a more compact region in the complex-2; plane 
than the zeros in Fig. [3[ We believe the zeros become dense in this region. Notice that this 
arch-shaped region is broader than the corresponding region in Fig. |2] because the turning 
points do not all lie along the same polar angle; therefore, the scaling fixes the magnitudes 
but not the positions of the turning points. Potentials with various values of a and 6 were 
also investigated and similar results were obtained. 

Example 3: ix 3 potential We obtain an ix 3 potential when we set N = 3 in (|1 . 1| ) . Using 
Runge-Kutta techniques in the complex plane, we have plotted the level curves of the real 
and imaginary parts of the complex eigenfunctions. By finding the intersections of these 
level curves, we have determined the zeros of the eigenfunctions numerically. These zeros, 
which are shown in Fig. |5|, lie along the Stokes' lines of the WKB approximation. Again, 
the zeros exhibit the complex version of interlacing. They have an imaginary part that 
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decreases as k increases. In Fig. |6] the zeros are scaled by z = x/\xtp\ = xE~ l l 3 , which 
fixes the magnitudes and positions of the turning points. This plot suggests that after the 
scaling the zeros become dense in the complex- z plane. Once again, this plot suggests that 
the distribution of zeros in the complex plane is a universal property of complex Sturm- 
Liouville eigenvalue problems associated with PT-symmetric Hamiltonians. 

Statement of the Conjecture 

From our studies we observe that the unsealed zeros of the complex eigenf unctions do 
not become dense on a contour or in a narrow region of the complex- x plane. In particular, 
for the ix 3 potential WKB theory predicts that ~ C/c 6//5 (k — > oo), where C is a constant. 
Thus, the turning points behave like xtp ~ (— iC) 1//3 /c 2//5 (k — > oo) and 

i /3r3/5 (2.2) 

ak 5 

Using Richardson extrapolation we have verified that the distance between zeros along the 
imaginary axis exhibits this /c-dependence. Consequently, the distance from the contour 
along which the zeros of lie to the contour along which the zeros of for large k lie is 
given by Y,T=o k~ 3 ^ which is infinitely far away. 

We have scaled the zeros by fixing the magnitudes of the turning points relative to a 
unit length. After this scaling is performed, the zeros appear to become dense in a narrow 
arch-shaped region in the scaled complex plane. If the zeros do become dense in this narrow 
region and exhibit the shifted complex version of interlacing, we conjecture that this behavior 
suggests that the eigenfunctions are complete in the scaled complex plane. Since we do not 
know what space within which to define completeness, we are unable to give a rigorous 
proof. 
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FIG. 1. Zeros of the 14th and 15th eigenfunctions, ^u(x) and ^i5(x), of the — (ix) N potential 
in the large- N limit. The plot shows the z-plane where the turning points have all been scaled to 
— 1 and 1 as in (2.1). The zeros lie in a small arch-shaped region in the complex plane. For any 
two zeros of Vl/15 (z), a zero of ^u(z) lies between them along the arch-shaped region. This is a 
complex version of interlacing. 
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FIGURE 2 
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FIG. 2. Zeros for the first 15 eigenfunctions of the —{ix) N potential in the large- N limit. The 
plot shows the z-plane where the turning points have all been scaled to —1 and 1 as in (2.1). The 
complex version of interlacing in Fig. 1 is again evident in this plot, but this plot also suggests that 
the zeros are becoming dense in a small region in the complex-z plane. 
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FIGURE 3 



-0.05 - 



-0.1 



-0.15 



-0.2 



-0.25 - 



-0.3 



-0.35 



-0.4 



□ □ 

■ ■ ■ 

o o 
• • • 

A A 



♦ ♦ ♦ 



o 

♦ o 



o 



» , 

t 

e s 

9 

e 

i_ 



♦ 

o * O 

, O O 



O « O 

O O 



O O O C5 o 

o o o o 



o 

o o 

O - ° 



off© 



♦ o 

. ° ♦ 



o » ° ° o e 
o » o o « 
o » „ ^ 



o a 
e" 
e 



zeros of Psi_1 (x 
zeros of Psi_2(x 
zeros of Psi_3(x 
zeros of Psi_4(x 
zeros of Psi_5(x 
zeros of Psi_6(x 
zeros of Psi_7(x 
zeros of Psi_8(x 
zeros of Psi_9(x 
zeros of Psi_10(x 
zeros of Psi_1 1 (x 
zeros of Psi_12(x 
zeros of Psi_13(x 
zeros of Psi_14(x 
zeros of Psi_15(x 
zeros of Psi_16(x 
zeros of Psi_1 7(x 
zeros of Psi_18(x 
zeros of Psi_19(x 
zeros of Psi_20(x 
turning points 



-1.5 



-0.5 



0.5 



1.5 



FIG. 3. Zeros for the 21 exactly solvable eigenfunctions of the quasi-exactly-solvable — x 4 po- 
tential in (1.2). The zeros lie in a small arch-shaped region in the complex-x plane and exhibit 
a complex version of interlacing. The zeros do not become dense in a region of the complex-x 
plane because the zeros lie on the Stokes' line of the wave function (the curve along which the 
wave function is oscillatory). These curves are different for each eigenfunction; the curves move 
downward as the energy increases and there is no limiting curve. 
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FIGURE 4 
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FIG. 4. Zeros for the 21 exactly solvable eigenfunctions of the quasi-exactly-solvable — x 4 po- 
tential (1.2) with the magnitudes of the turning points fixed. The plot shows the complex-z plane 
in which the magnitudes of the turning points are fixed to unit length. In the z plane the zeros 
lie in a much narrower arch-shaped region than in Fig. 3. Once again, the zeros exhibit a complex 
version of interlacing, and now, we believe they become dense in a narrow region of the scaled 
complex plane. 
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FIGURE 5 
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FIG. 5. Zeros of the first six eigenfunctions of the ix 3 potential. These zeros lie in a small 
arch-shaped region in the complex plane and exhibit a complex version of interlacing. The zeros 
do not become dense in a region of the complex plane because, as in Fig. 3, the zeros lie on a 
sequence of curves that move downward with increasing energy and remain well separated. 
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FIGURE 6 
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FIG. 6. Zeros for the first six eigenfunctions of the ix 3 potential with the turning points fixed. 
The plot is done in the complex- z plane where the magnitudes of the turning points are fixed to 
unit length. As a result, the zeros lie in a narrower arch-shaped region in the complex plane than 
in Fig. 5. The zeros exhibit a complex version of interlacing and we believe they become dense in 
the scaled complex- z plane. 
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